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ABSTRACT 

Context. Despite all the studies, the geometry of the wind at the origin of the blueshifted broad absorption lines (BAL) observed in 
nearly 20% of quasars still remains a matter of debate. 

Aims. We want to see if a two-component polar+equatorial wind geometry can reproduce the typical BAL profiles observed in these 
objects. 

Methods. We built a Monte Carlo radiative transfer code (called MCRT) to simulate the line profiles formed in a polar-l-equatorial 
wind in which the photons, emitted from a spherically symmetric core are resonantly scattered. Our goal is to reproduce typical C iv 
line profiles observed in BAL quasars and to identify the parameters governing the line profiles. 

Results. The two-component wind model appears to be efficient in reproducing the BAL profiles from the P Cygni-type profiles to the 
more complex ones. Some profiles can also be reproduced with a pole-on view. Our simulations provide evidence of a high-velocity 
rotation of the wind around the polar axis in BAL quasars with non P Cygni-type line profiles. 

Key words. Quasars: absorption lines - radiative transfer - Methods: numerical 



■ 1. Introduction 

Depending on the selection technique and the definition used, 
about 20% to 30% of the quasars detected in recent surveys show 
the presence of the broad absorption line (BAL) troughs associ- 
ated with the emission lines in their rest frame UV spectrum (e.g. 
Knigge et al. 2008, Ganguly et al. l2008l l. These BALs, reminis- 
cent of the P Cygni-type profiles seen in the spectra of massive 
stars, are mainly observed in high ionization lines like C iv and 
Si IV and are sometimes detected in lower ionization species like 
Mg II. They reveal strong outflows from quasars (Scargle [19721 ). 
which can reach velocities up to 0.2 c (Foltz et al. |1983l ). 

Despite the large number of observations, the physical and 
'] geometrical properties of the wind at the origin of the BALs 
. remain largely unknown (e.g. Brotherton I2007I) . Moreover, the 
' distance at which those objects are found (z > 1 .5, so that C iv is 
' shifted in the optical domain) hampers direct observation of the 
regions at the origin of the BALs even with the best telescopes 
presently available. Thus, all the information we can get about 
, the inner regions of BAL quasars comes mainly from indirect 

■ observations. 

The first attempts to model the BAL profiles considered the 
resonant scattering of photons emitted by a continuum source 
in a spherically symmetric stellar-like wind (e.g. Scargle et al. 
119721 Surdej & Hutsemekers I 19871 ). However, the growing num- 
ber of observed spectra displaying a huge variety of line profiles 
(Korista et al. 119931 ) revealed the need for other wind models. 
Facing the diversity of line profiles, Turnshek ( |1984al l proposed 
that BAL quasars could be broadly divided into two samples: 
those quasars that exhibit smooth P Cygni-type profiles, and 
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Other ones that display an absorption trough that is detached in 
velocity from the associated weaker and wider emission peak. 
These observations indicate that the properties of the wind are 
more complex than the simple spherically symmetric outflow 
inferred for stellar winds (Lee & Blandford |1997l l. However, as 
emphasized by Turnshek (1984b ), it is very likely that distinct 
types of BAL QSOs do not exist but are instead different mani- 
festations of the same phenomenon. 

The similarities of the emission line, optical continuum, and 
infrared properties of BAL and non-BAL QSOs (e.g. Weymann 
et al. 1991 , Gaflagher et al. [T999l Reichai'd et al. 2003 Gaflagher 
et al. 2007 ), as well as the spectropolarimetric observations (e.g. 
Schmidt & Hines 1999 Ogle et aL .1999 , Lamy & Hutsemekers 
2004), favor a unification by orientation scheme for the BAL 
QSOs over the evolutionary scheme (Hazard et al. 1 19841 Becker 
et al. 1 20001 ). In the unification by orientation scheme, only a frac- 
tion (roughly corresponding to the observed fraction of BAL 
QSOs) of the continuum source is covered by optically thick 
material producing the broad absorption lines, which suggests a 
disk-like equatorial geometry for the BAL region (e.g. Turnshek 
1984a Hamann et al. 1993 Murray et al. [19951 Elvis et al. 
2000 Yamamoto 2002). Such a geometry is supported by the- 
oretical studies and commonly accepted, since the QSOs are 
thought to be powered by accretion of matter onto a supermas- 
sive black hole in the form of a disk, from which the wind could 
be launched. However, the recent discovery of radio loud BAL 
QSOs (e.g. Becker et al. 120001 ) and subsequent radio variability 
studies reveal polar outflows in at least some of them (Brotherton 
et al. 12006) Zhou et al. 12006) Ghosh & Punsly |20D71 ). Models 
combining polar and equatorial components have also been sug- 
gested (e.g. Lamy & Hutsemekers 2004) and evaluated from a 
theoretical point of view (Pereyra et al. 120041 Proga et al. 120001 
Proga|2003i Proga & Kallman|2004l) 
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In this context and given the similarities between typical 
BAL profiles (e.g. Korista et al. '1993) and the line profiles pro- 
duced by a two component polar+equatorial wind like the one 
presented by Bjorkman et al. (1994), our goal in this first pa- 
per is to determine whether such a simple two-component wind 
can qualitatively reproduce the various types of line profiles ob- 
served among the BAL QSOs. We also try to identify the key in- 
gredients needed to reproduce BAL profiles. In a second paper, 
we will investigate the effect of microlensing on these profiles, 
aiming at a realistic interpretation of the spectral differences ob- 
served in gravitationally lensed BAL QSOs like H1413-H1 17 (cf. 
Hutsemekers et al. 2009). 

In Sect|2l we present MCRT, the Monte Carlo radiative trans- 
fer code we implemented in order to simulate resonance line 
profiles in a two-component axisymmetric wind. In Sect|3] we 
briefly identify the influence of the wind model parameters on 
the line profiles we computed. In SectlH we show how MCRT 
is able to reproduce typical C iv BAL QSOs line profiles. We 
discuss the results of the line profile fitting and summarize our 
conclusions in the last two sections of the paper 

2. The MCRT code 

MCRT is a Fortran77 fully 3D Monte Carlo (MC) radiative 
transfer (RT) code that we built to compute the resonance line 
profiles produced in axisymmetric winds. The use of the Monte 
Carlo simulation technique allows the radiative transfer equation 
to be solved exactly (i.e. without making use of the Sobolev ap- 
proximation), as well as ensuring the self consistent treatment of 
the radiative coupling between distant regions in a wind subject 
to more complex velocity fields than monotonic radial laws (e.g. 
Kniggeet al. 1995). 

Monte Carlo RT code have been extensively described (e.g. 
Knigge et al.[l995, Wood et al. l200n Dijkstra et al. 2006), so that 
we only recall here the fundamental principles of this technique 
and the particularities of the code we developed. 

As stated in the introduction, our main goal is to identify 
of the key ingredients (geometry and overall kinematics) of the 
wind governing the typical profile of the BAL QSO UV reso- 
nance lines. Thus we do not consider negligible effects, such 
as the relativistic ones that remain small even for outflows with 
high (Vmax ^ 0.2c) terminal speed (Hutsemekers & Surdei [T990] l 
or the fact that the line is a resonance doublet, because the veloc- 
ity separation of the doublet components is small with respect to 
Vmax (e.g. Hewitt et al. |T974l Grinin |T984l ). 

2.1. Radiative transfer with Monte Cario tecliniques 

When using the MC technique, the solution of the RT equa- 
tion is found by following a huge number of photons on their 
way through the wind. Each step in the photon's life (position 
and direction of emission, position of interaction, etc) is deter- 
mined by the mean of random numbers distributed according to 
the normalized probability density function (NPDF) of the cor- 
responding simulated physical process. Thus if one wishes that 
the frequency v, of all of the emitted photons follows a given 
law L(y) over the frequency interval {Vmin,Vmax\, then v, will be 
randomly chosen by solving the transformation equation (Press 

et ai. fiyga : 



I 



L{v)dv, 



(1) 



where ^ is a random number drawn from a uniform distribution 
in the interval [0, 1 ] . In MCRT this number is generated using the 



"ran2" subroutine of Press et al. (I1992I I. In the following, each 
new occuiTence of ^ refers to the call of such a new random num- 
ber. In general, there is seldom an analytical solution to EqUlso 
we implemented the "table lookup method" (see Avery & House 
1968), which allows us arbitrary NPDF's. 

In MCRT the initial position of emission of the photons is 
chosen isotropically on the surface of the continuum emission 
region, which is modeled by a sphere of radius and of infi- 
nite optical depth {tc - oo) located at the center of the wind. 
The direction of travel through the wind is then determined by 
randomly sampling a half sphere taking into account that the 
photons are forced to leave the continuum source upward. 

2.2. Continuum ptiotons and resonance scattering 

The wind is filled with 2-level atoms whose rest-frame normal- 
ized absorption profile (pahs is described by a Gaussian (Natta & 
Beckwith 1986. Knigge et al. il995i) such that 



PahAv - Vo) = 



(pabAv - Vo, CTturb) if |v - Vol < \^Vah,\ 

everywhere else 



(2) 



where vo is the rest-frame frequency of the considered transition, 
and 



<t>abs{v- V0,O-,urb) 



1 



exp 



urb 



-(v - Vo) 



2^ 



lurb 



+ K\ , (3) 



in which K ensures a continuous transition between the absorp- 
tion profile and the zero intensity. Here, // is a constant allow- 
ing for the normalization of the line profile over the interval 
Vo ± |Avaij|/2, where Ava^j is the full width at zero intensity 
(FWZI) of the absorption profile. We choose the value of hvabs to 
ensure the continuity of the absorption profile at the border of the 
interval vo+|Avais|/2. The parameter cr,„ri - Av,„ri/(2 V2 In 2) is 
such that IS.v,urb = 2 vo(vtui-b/c) is the FWHM of the absorption 
profile. The velocity Vmrb includes the thermal and the macro- 
scopic turbulence components in the wind that broaden the ab- 
sorption profile. We assume Vtmb to be constant throughout the 
wind. 

Owing to the velocity field "\^(r, 6, (p) present in the wind, the 
initial frequency v,- of a photon flying in the direction it is seen 
Doppler-shifted by an atom of the wind in such a way that its 
"local" frequency v; in the atom rest-frame is given by 



V/ 



^1? 



(4) 



Thus a photon will enter in resonance with the surrounding 
atoms only if its local frequency fulfills the condition defining 
the so-called "resonance zone": 



AVfli.5 AVabs 
Vo r— < V, < Vo H- — -— 



(5) 



When the photon enters such a region, the opacity of the medium 
becomes nonzero as does the probability of being absorbed. If n 
resonance zones are found along the direction of propagation of 
the photon in the wind, the total optical depth r,o, seen by the 
photon until it escapes the wind is simply computed as 



TtotiVi) = Kva4>abs{vi ' Vo, CTturb) ds , 

7=1 



(6) 



where Ky^ is the total absorption coefficient of the considered res- 
onance transition, and aj and bj are respectively the coordinates 
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of the beginning and of the end of the / resonance zone found 
along the line of flight of the photon. 

Given the probabilistic interpretation of the RT, a photon 
experiencing a total optical depth of t,o,(v,) has a probability 
p - g-T,„,(vi) Qf escaping the medium without being absorbed. 
This interpretation is used in the MC code to identify the occur- 
rence and the position of the scattering sites along the path of 
the photon through the wind. Indeed by using the transformation 
equation (Eq[TJ we can determine the random optical depth tmc 
at which the photon interacts : 

TMc = -ln(l-^). (7) 

Because we stored the run of r,of = Tfof(s) along the photon 
path, it is easy to invert this relation and find the point s where 
Tto/ - Tmc- At this location the photon is radiatively absorbed 
and then instantaneously re-emitted at a frequency and in a direc- 
tion chosen by assuming a complete redistribution in frequency 
and direction (CRFD, Lucy 1971 , Mihalas et al, 1976). This pho- 
ton may then either be re-absorbed somewhere else in the wind 
or leave it and be detected by one of the detector (spectrographs, 
imagers) located around the wind. 

To decrease the simulation time in the case of non- 
spheric ally symmetric winds we make an intense use of the ad- 
vanced concepts of "first forced interaction" (e.g. Cashwell & 
Ev erett[T959[ Witt jlWTI l and "peeUng off" (e.g. Yusef-Zadeh et 
al. II984I Wood & Reynolds 11999) where we follow a photon 
packet rather than a single photon. 

We checked the validity of our MCRT code by compar- 
ing the line profiles we obtained to line profiles computed with 
two traditional methods for spherical winds that allow an ex- 
act integration of the transfer equation. These benchmarks are 
the well-known SEI method of Lamers et al. (II987I I and the 
comoving frame method of Hamann et al. (1981). We noted 
good agreement between the general shape of the computed pro- 
files, regardless of the considered turbulence (AF/F < 5% on 
the normalized emission peak flux, as well as a good match 
between the absorption profiles). We also tested MCRT in the 
case of axi-symmetric winds by comparing the profiles obtained 
with MCRT to those produced by the SEI method adapted by 
Bjorkman et al. (1994). Once again we observed good agree- 
ment between the line profiles produced by both methods (see 
Borguet 2009 for details). 

2.3. Pure emission 

While the shape of the C iv line in BAL QSOs can be mostly gov- 
erned by resonance scattering (Scargle et al. 1 19721) . the presence 
of C III] emission constitutes evidence that part of the emission 
is due to collisional excitation (Turnshekl 1 984a Turnshek il988i 
Hamann et al. ll993l l. To account for this second source of pho- 
tons, we allow the production directly in the wind of a fraction of 
photons fe = I pyre emission I Iconiinuum- The choicc of the location of 
the emission (r^. Be, (f>e) of these photons is made using a random 
sampling of the corresponding NPDF: 

H>e rO, nr, 

^ ^ p(re,0e,<f>e) ^ \ I I T](r, 6, (f>)r^ siu 6 drdOd^) , (8) 
Jo Jo Jr;„ 

where the ri{r,0,^) is a function that describes the emissivity 
throughout the wind. Once again, our goal here is not to provide 
a detailed self-consistent model of the wind so we choose as a 
first guess an emissivity function of the form 

IR V 

T]{r, 0, 0) = n(r, 6, 0) , (9) 



where n(r, 0, cp) is the density of the ion through the wind and 
where the second term allows taking the temperature distribu- 
tion and the ionization fraction into account. In the following we 
simply take y = 1 in order to reduce the number of free parame- 
ters. 

Each photon emitted then makes its way through the wind 
where it can be scattered and then finally escapes the wind in 
order to be detected by a distant observer. 



2.4. Tine wind modei 
2.4.1. Introduction 

As stated in the introduction, it is diflScult to give a simple ex- 
planation to the observed BAL profiles when using a spherically 
symmetric expanding wind. The obvious next type of geome- 
try that can then be considered is the axi-symmetric one. The 
simplest of these models would still consist in a wind origi- 
nating from the central core. Such a generic model, with polar 
and equatorial components like the one presented in details by 
Bjorkman et al. (119941) produces line profiles remarkably simi- 
lar to those observed in some BAL QSOs. Although simple, this 
model is versatile enough to produce a variety of line profiles, 
as observed in BAL QSOs. It constitutes a first good approxima- 
tion to the more complex wind from disk models proposed for 
AGN outflows in which the BALR and the BELR are generally 
cospatial (Sect.[T]i. 

In our model, we adopt stellar wind laws to describe the 
kinematics of the winds observed in quasars. Indeed, quasar 
winds are also supposed to be driven by radiation (Arav & Li 
1996 Murray et al. 119951 ) as suggested by the line-locking in 
the spectra of some BAL QSOs (e.g. Weymann et al. 119911 
Korista et al. 119931 Arav 119961 II997I ). However, there are im - 
portant diff'erences between stars and quasars (e.g. Aravl l994] l. 
One of them is the so-called overionization problem caused by 
the strong UV/X-ray central source in quasars (e.g. Proga et al. 
2000). Several scenarios have been suggested to solve this prob- 
lem: Mun-ay et al. (1995), Murray & Chiang (1997), Risaliti & 
Elvis (2009), Punsly ( 1999) and Ghosh & Punsly (120071) . In our 
study we assume the existence of shielding material between the 
radiation source and the outflow that prevents the total ionization 
of the outflow (see Krolik ll999b . 

Another difference between stellar objects and quasars 
comes from a significant fraction of the radiation in quasars sup- 
posedly being emitted from an accretion disk rather than from 
a spherically symmetric photosphere (e.g. Proga et al. 20001 1. 
Outflows with axial geometries have been studied by several au- 
thors, who show that the flow can be launched vertically from 
the disk and then pushed away by the radiation from the central 
source (Murray et al. 119951 Proga et al. I2000i ) with the eleva- 
tion of the wind over the disk is still small when the flow starts 
to expand radially. Since the launching conditions are unclear 
(Risaliti & Elvis 2009) and since we assume the BALR and the 
BELR to be cospatial, we simplify the geometry by assuming 
that the wind is purely radial, which is correct at some distance 
from the disk. We considered a wind launched from a sphere 
with a radius equal to the inner radius of a typical BELR (10^^ 
pc). The outer radius of the BELR/BALR was chosen to be 1 
pc, the radius at which the wind reaches the terminal velocity. If 
needed, the anisotropy of the radiation field from the continuum 
source can be readily introduced into our model. 



4 



B. Borguet & D. Hutsemekers: Fitting the C iv BAL profiles in quasars with a poIar+ equatorial wind model 



z ^ 



/ ^ i 



/ I'lim 



\ 

! A9 

_ _ Observer 

Fig. 1. Illustration in the x - z plane of the two-component wind 
used in our MC simulations, the velocity field is pictured by 
arrows in the lower part of the picture. The equatorial wind is 
shown edge-on (in dark grey) so that the wind model is rotation- 
naly symmetric along the z axis. The viewing angle of the ob- 
server is / = 90°, and the absorption profile is produced by both 
component, the polar one contributing only at low velocities (i.e. 
for X < riiin) given the disk opening angle and the viewing angle. 
Inspired from Fig. 8 of Bjorkman et al. 019941 1. 



2.4.2. Velocity field and density law 

Keeping in mind the aforementioned simplifications we then 
decided to implement the model described by Bjorkman et al. 
d 19941 1. We recall its basic characteristics here. The velocity field 
'^{r,0,4>) considered here can be written in spherical coordi- 
nates: 



Vr(r, 0)'e*r + wAr, 0)ei,+ Ve(r, ff) el , 



(10) 



where for simplicity we assume, as in Bjorkman et al. (1994), 
that the streamlines lie on surfaces of constant polar angle so 
that \q - 0. The radial component of the velocity field has the 
typical jS-law shape: 



Vr(r, 6) = Vmin "I- (Vniax(6') - Vmin) | 1 - 



(11) 



where v^in is the wind speed at the surface of the source of the 
continuum. The variable terminal speed Vm-ixiQ) allows for the 
introduction of a slow expanding equatorial wind in a faster polar 
wind: 



x(0) 



+ (V„ 



(12) 



where v^^^ ™d vJJ^ax respectively the terminal velocity of the 
equatorial and the polar components. The function that allows 
the smooth transition between the two components of the wind 
takes the form 



m = 0.5 



1 + - arctan 



sinAO - cos 6 



A0, 



(13) 



where AO is the equatorial wind half opening angle and where the 
parameter A0, controls the transition between the equatorial and 
the polar component of the wind. We chose, as in to Bjorkman 



et al. (119941) . a low value for A6, = 0.001 in order to keep this 
region much thinner than the equatorial wind opening angle. 

The v^ component of the velocity field is simply given by 
assuming conservation of the angular momentum in the wind: 



Vfot Rin . „ 

V0 = sm 6 , 



(14) 



where Vrot is the rotational speed at the surface of the source of 
the continuum. 

The law governing the distribution of the ion density is de- 
rived from the equation of continuity of matter in a central wind 
and is parameterized using 



\t<i„ 



Vr(r,0) 



(15) 



with the parameter a implicitly including the radial variation of 
the ionization (i.e. a = 2 meaning a constant ionization through- 
out the wind) and «()(0) allowing the transition between the polar 
and the equatorial components: 



«o(0)=nr+(«o'-«r)/W' 



(16) 



where n^^ and are, respectively, the density of the considered 
ion at the base of the wind in the equatorial or in the polar com- 
ponent. The value of n^" is computed by specifying the value 
Ttotp of the total optical depth of the wind along the polar axis 
(i.e. 9 = 0°) integrated over frequency: 



XAv, 



T,o,(Vi)dVi , 



(17) 



with 2 Av/ the measured width of the observed line profile. The 
value of is fixed by the free parameter kp,,,, which defines 
the ratio of the ionic density between the two components at the 
base of the wind: 



po 



(18) 



3. Parameter study 



Here we concentrate on the main parameters affecting the line 
profiles for a two-component wind where we suppose no pure 
emission (i.e. fe = 0). We do not discuss the effects of the pa- 
rameters governing the velocity and density laws or the effects 
of the turbulent component in the wind since they are similar 
to those observed in the well-understood spherically expand- 
ing wind case (Castor & La mers[T979 , Beckwith & Natta ll987l 
Hamman 1981 Lamers et al. |1987l l. 

The most important parameters we have to specify for com- 
puting a line profile are the frequency integrated polar optical 
depth T,o,^_, the ratio between the equatorial and the polar ionic 
density kp„, , the velocity ratio between the polar and the equa- 
torial terminal speed Vm-dx/^max^ the disk half-opening angle AO, 
the viewing angle /, and the ratio 'VroJ^mux between the rota- 
tional speed of the source of continuum and the polar termi- 
nal velocity. The parameters used for the reference line profile 
are summarized in Table [T] Such a parameter study was previ- 
ously carried out by Bjorkman et al. (1994) using an SEI-type 
method. Because their calculations agree with ours using an ex- 
act method, we only recall here the general effects produced by 
each parameter and refer the reader to the Bjorkman et al. pa- 
per for further details. We do, however, emphasize the effect of 
the wind rotation given the important changes in the line profile 
produced by the variation of this parameter, more particularly 



B. Borguet & D. Hutsemekers: Fitting the Civ BAL profiles in quasars with a polar+ equatorial wind model 



2.0 



1.0 



1 1 I 1 




lilt 





-0 

1.0 



0.0 






■ 1 I 1 ■ 






■ ki, = 35tKHiB^ 




■ 1-2:,. »5(incn™,'> 













-2.2 -0,8 -OA 



0,0 

X 



0.4 




0.3 1,2 



.a 1.0 



0.0. 
2.0 



I.O 



1 1 1 1 1 1 . 


( ' r- ' 






i = 7J- 




- i = «ir I ' 




l=S?' / 




1 









0.0 















' • iilXU-wAi 












V„ = TJ(ICIi-jn/i 








1. 





-1.2 



0,0 
X 



OA 



as 



Fig. 2. Illustration of the effects of the wind's main parameter on the line profile (see text for details). Each spectrum is represented 
as a function of the normalized wavelength X - (A- A^))/{A„,„^ - Aq) and is the result of the simulation of 2 10^ photon paths through 
the wind (S /N > 100 on the continuum level). 



Table 1. Wind parameter for the benchmark model 



Parameter Value 


Parameter Value 


^ max 


10000 km s"' 


a 


2.00 


^min 


100 km s"' 


P 


1.00 


'max 


2500 km s"' 


Ae 


15° 


Vrot 


Okms"' 


AO, 


0.001 


Vturb 


500 km s"' 


i 


90° 




4.00 


fe 


0.0 




10.0 


J 


1.0 



when the equatorial disk is viewed near edge-on (Mazzali |1990l 
Petrenz & Puis 1996 Busche & HillierEoSj- 

The profiles illustrated in the panel A of Fig|2]are represen- 
tative of the profiles that can be produced in a two-component 
wind when the equatorial component is seen near edge-on. The 
shape of these profiles is constituted by typical P Cygni-type 
profile extending to the higher velocities coming from the polar 



component, but where a sharp absorption trough is produced at 
lower velocities by the slowly expanding equatorial wind. In this 
panel we illustrate the evolution of the line profiles as a function 
of the total polar optical depth integrated over frequencies Tiot,,- 
Similar to the spherically symmetric expanding wind case (e.g. 
the atlas constructed by Castor & Lamers 1979), when T,ot^ is 
increased, we observe an increase in the emission peak and in 
the equivalent width of both the emission and absorption parts 
of the line profile. No variations are observed in the equatorial 
absorption since that part of the profile is already saturated for 
the lower value of t,o,^. 

The effect of the kp,„ ratio is illustrated in the panel B of 
Fig. |2] When kp,„ is decreased, the depth of the absorption pro- 
file of the equatorial component evolves in the same way. As 
in Bjorkman et al. (1994), we observe that the equatorial com- 
ponent is still optically thicker to the radiation than the polar 
component even for kp^ =0.1, and this because of the narrower 
velocity range of the equatorial component. 
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The change in the ratio between the teminal speed v^ax ™d 
Vmax of the two Components essentially modifies the position of 
the edge of the low-velocity absorption component (see panel 
C of Fig. |2|. We also observe that, when v^j^^x is increased, the 
absorption due to the presence of the disk no longer remains 
black near the line center. Indeed, to produce a black absorption, 
the disk has to cover the whole source of continuum. But for an 
equatorial wind with a half-opening angle A0, this only happens 
beyond the radial coordinate r > r/„„ = cotanAfl (see Fig.[Tl), 
which corresponds to the velocity vjj^ ~ Vr(riin,, 90°), the value 
of which is sensitive to \m-dx{d) (cf. Eq.fTTI). 

In panel D of Fig. [2] we illustrate the evolution of the line 
shape as a function of the half-opening angle of the equatorial 
wind A6'. We observe the development of a broad high-velocity 
absorption when A9 is decreased. This observation is once more 
explained by considering the radial coordinate r/,„, at which the 
disk completely covers the source of the continuum. As the high- 
velocity absorption (due to the polar component) occurs before 
the equatorial wind, the high absorption velocity range is limited 
within the interval [0, vj'j^], where v^"^ ~ Vr(rii,n, 0°). In contrast, 
when AO is increased, we observe the apparition of a secondary 
blueshifted emission peak beyond v^^^^ given that, for a suffi- 
ciently wide disk, there is no remaining polar absorption of the 
continuum. 

The viewing angle / to the wind plays a critical role in the 
line profile (see panel E of Fig. |2]i. Indeed the equatorial ab- 
sorption component is only seen when the disk is viewed near 
edge-on (i.e. / = 90°). Because / is decreased from an edge-on to 
a pole-on (/ = 0°) view of the disk, the low-velocity absorption 
produced by the disk decreases, since only the high- velocity end 
of this component covers the continuum. 

Finally, panel F of Fig. |2] illustrates the modifications of 
the line profile when varying the rotational speed of the wind. 
The major effect of the rotation is observed on the emission 
peak, which decreases and is significantly redshifted when the 
V,ot/Vmax ratio is increased. Indeed, when there is no rotation, 
the emission peak is relatively sharp because that part of the pro- 
file is produced in the inner, optically thicker regions of the wind, 
whose velocity range is narrow. When rotation is considered, the 
resonance zones become twisted in such a way that the inner re- 
gions are distributed over a wider range of projected velocity. 
This induces the smoothening of the emission peak since the ro- 
tation allows the absorption to occur in the redder part of the 
continuum at frequencies centered on the rest-frame frequency 
of the line transition vo (e.g. X - 0.0) (see Figl3]and also Hall 
et al. [20021 for a detailed explanation). The bluer part of the ab- 
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sorption profile is not significantly affected by the rotation be- 
cause this part of the profile is produced in front of the contin- 
uum source, i.e. where the projected rotational velocity is nearly 
zero. 
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Fig. 3. Closer illustration of the effect of the rotation of the wind 
over the line profile, its part in emission, and its part in absorp- 
tion. The dotted line represents the profiles for which the rota- 
tional velocity is zero, while the straight line represents the pro- 
files with Viot - 0.25 v|^nax' ^nd the dashed line a wind where 
Vrot - 0.50 Vmax- The parameter used are those gathered in 
Table[T](see text). 



We retained as far as possible objects free of doublet effects 
(i.e. » Civ doublet separation), non-blended and with 
smooth emission/absorption C iv line profiles (i.e those less 
affected by multiple absorption troughs). 
We selected objects for which the continuum can be securely 
defined on both sides of the C iv line profile, allowing for a 
correct normalization of the spectra. 



4. Fitting BAL profiles with the polar+equatorial 
wind model 

In this section we try to reproduce a representative sample of 
C IV line profiles observed in BAL QSOs and extracted from the 
paper of Korista et al. ( 1993), which provides homogeneous and 
high-quality spectroscopic data for a large sample of BAL QSOs. 
Our goal is to check that such profiles can be modeled using a 
simple two-component wind geometry for the BALR. Ten spec- 
tra were carefully selected out of the 72 BAL QSOs presented in 
the Korista et al. paper on the basis of two criteria. 



The observed and the fitted spectra of each selected object 
are illustrated, normalized to the continuum level, in Figs.|4]and 
|5] We note good qualitative agreement between the observed line 
profile and the model. The parameters used in the fit of each 
C IV line profile are collected in Table |2] While some parame- 
ters (v|j^ax' ^mux^ U '''toi,,, /«■■■) wcrc allowcd to vary, we tried to 
keep fixed the radial velocity gradient in the wind (the param- 
eter), the emission parameter y and the disk opening angle AO. 
The latter was chosen equal to the value generally adopted for 
the opening angles of disk winds in AGN models (Murray et al. 
1995, Proga & Kallman 2004). However, even those parameters 
were varied in a few cases to slightly modify the part of the pro- 
file in emission. In particular, we found that the opening angle 
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M of the disk has to be decreased for fitting the high-velocity ab- 
sorption component (see Sect. |3]l observed in the P Cygni-type 
Une profiles (Q1333-H2840 and Q1413-H1143). 

The sample contains various types of BAL spectra, from the 
P Cygni-type to the more complex profiles discussed in the paper 
of Korista et al. ( I1993I I. This subdivision of BAL QSO spectra 
was introduced by Turnshek (1984b) and based on the visual 
properties of the BAL profiles, where the P Cygni-type profiles 
with a smooth absorption and a high emission peak relative to the 
continuum level are distinguished from the complex/detached 
absorption ones with a lower emission peak. However we tried 
to model all these line profiles using the same wind geometry 
since, as suggested by Turnshek (, 1984b ). it is likely that distinct 
types of BAL profiles are the consequence of a unique mass-loss 
phenomenon. 

The empirical fitting procedure we adopted here makes the 
distinction between three types of line profiles, each possessing 
its own set of spectral signatures. The first kind of line profile is 
the P Cygni-type (Q1333+2840 and Q1413+117 see top pan- 
els of FiglH). In this subsample, the absorption trough seems 
constituted of two subtroughs covering two overlapping veloc- 
ity ranges, i.e. an optically thick "narrow" (~ 4000 km s"' wide) 
absorption component superimposed on an optically thin absorp- 
tion trough extending up to velocities of 10000 km s"'. These 
characteristics suggest that there is an equatorial denser disk-like 
region seen nearly edge-on and producing the narrow absorption 
observed. The necessity of considering such a more slowly ex- 
panding disk in a polar wind is illustrated particularly well when 
observing the C iv line profile of Q1413+1 17. In this object, the 
absorption profile extends to velocities up to -10000 km s"', 
which implies, through resonance scattering, the presence of an 
underlying emission component extending from -10000 km s ' 
to 10000 km s '. Since the flux is almost zero at slower veloc- 
ities (v > -4000 km s"'), a part of these photons has to be re- 
absorbed somewhere, suggesting there is an optically thicker re- 
gion on the line of sight to the observer This observation fits the 
framework in which a part of the emission from the BELR is 
absorbed in the BALR (e.g. Turnshek et al. |I988l Hamann et al. 
[1993). 

The second type of line profile we considered is presented 
in the lower panels of Fig. |4l These profiles are character- 
ized by an asymmetric emission peak whose intensity rela- 
tive to the continuum level is significantly lower than in the 
P Cygni-type ones. Some individuals of this subsample also 
show, similar to the P Cygni-type ones, higher velocity ab- 
sorption through superimposed to a narrower component (best 
seen in, e.g., Q0842 -1-3431, Q1235-i-1453). The introduction of 
a significant rotation of the wind around the polar axis is re- 
quired to produce the broad asymmetric emission peaks ob- 
served in the C iv line profile of seven out of the ten objects in our 
sample (Q0019+0107, Q0145-I-416, Q0226-1024, Q032-3344, 
Q0842+3431, Q1235+1453, and Q1239-I-0995). The rotational 
speed generally needs to be quite higher than the polar termi- 
nal velocity (reaching up to 50 % of v^^^^) and is approximately 
equal to the equatorial wind terminal velocity. For the P Cygni- 
type profiles (Q1333-I-2840, Q1413+1 143 or even Q0041-4023), 
it is also necessary to consider a rotation of the wind around the 
polar axis, although it must be smaller to properly account for 
the higher emission peak and the sharp transition between the 
emission and the absorption parts of the profiles. 

In each case, the best fit was chosen by eye, once the model 
was qualitatively similar to the profile of the observed line. The 
simulation time needed to produce a single line profile in these 
optically thick winds prevents us from using a.x^ type technique 



while searching for the best model. However, this is not a main 
drawback since our main goal is to show that a simple wind 
model is able to approximately reproduce a variety of resonance 
line profiles observed in BAL QSOs. Moreover, given the degen- 
eracy between some of the model parameters {kp,„ with /, v^^ax 
with /, j6 with a, etc), as well as the difficulty in some cases of 
correctly evaluating vlJ^^x or other model parameters, more than 
one best fit is generally possible. 

The absence of the two-component wind signatures defined 
above in the C iv line profile of some objects led us to define a 
third subsample of line profiles. These line profiles can be fitted 
by a two-component wind (see upper panel of Fig. |5]). However, 
some profiles (the prototype being Q0019-1-0107) do not show 
evidence of a polar absorption component at high velocities, sug- 
gesting that the polar outflow may not be present in these objects. 
Several tests performed on such line profiles with MCRT showed 
that Q0019-i-0107-type line profiles can be produced in a single, 
rapidly rotating equatorial wind seen nearly edge-on (see lower 
left panel of Fig.jSj- In the same vein, in the bottom right panel of 
Fig.|5] we illustrate how the C iv line profile of Q0041-4023 can 
be reproduced by a two-component wind seen nearly pole-on. 
This arises when a single deep absorption trough is associated 
with a quasi-symmetric emission profile. Once again, given the 
uncertainties on the wind parameters because of the lack of clear 
signatures in the line shape, several wind models can be fitted to 
the observed line profiles (here we chose as typical values for 
model M2 v^^^^ = = 0.5 vl^"^). 

From Figs. |4]and|5] we note that our simple model is able 
to reproduce the diversity of the BAL profiles observed in a real 
sample of objects quite well. Interestingly enough, in order to be 
able to fit the C iv profiles with MCRT, we must shift the whole 
simulated line profile with respect to the emission peak, usually 
used for redshift determination. This shift is needed to center the 
underlying emission component of the profile on the zero veloc- 
ity. Indeed, when dealing with resonant scattering, an absorption 
trough extending over the velocity range [-v^ax, 0] produces an 
emission feature extending from -Vmax to -nvmax- The redshift 
of the quasar determined from the center of the underlying C iv 
emission line is given in the second column of Table|2] 

From the fitting procedure, one of the major parameters that 
control the line profile in this type of wind remains the viewing 
angle /, which plays a crucial role in the shape of the absorption 
part of the line by controlling the relative contribution of the 
equatorial and polar components (when both are required). Thus 
when a line profile exhibits a sharp, deep absorption trough su- 
perimposed on a shallower high-velocity absorption component, 
the quasar is probably viewed along a line of sight, such as the 
dense equatorial wind seen nearly edge-on (e.g. Q1413-I-1 17). 
However, it should be kept in mind that an edge-on disk does 
not necessarily produce a completely black absorption trough, 
even for the highest optical depths (cf. Table |2]i, since photons 
scattered to the observer may not be completely reabsorbed by 
the disk then filling in the absorption trough. This agrees with 
what is usually observed among the Korista et al. ( 1993) spectra 
and has already been pointed out by Lee & Blandford (il997j or 
Arav et al. (120071 and references therein). 

Finally in all modeled profiles, we need to allow for the cre- 
ation of photons inside the wind (a fraction fe > relative to the 
continuum intensity). This intrinsic emission produces stronger 
emission peaks than in the case of a purely resonant scattering 
wind. 
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Table 2. Parameter of the two-component wind model for the best fit of the C iv BALs in the selected quasar sample. 



BAL QSO Name 


z 






Vmax (kms ') 


Vmax (kmS ') 


V,-„, (kms-') 


a 


P 




in 


fe 


y 


Q1333+2840 


1.897 


3.00 


10.0 


13000 


4000 


1000 


3.0 


0.8 


10 


80 


0.21 


1.0 


Q1413+1143 


2.522 


10.0 


15.0 


11000 


3500 


2000 


3.0 


0.8 


10 


88 


0.16 


1.0 


Q0145+0416 


1.990 


1.00 


12.0 


18000 


9600 


6400 


2.5 


1.2 


15 


80 


0.18 


1.0 


Q0226-1024 


2.164 


1.00 


10.0 


25000 


12500 


12500 


2.0 


1.5 


15 


87 


0.10 


1.0 


Q032 1-3344 


1.955 


10.0 


10.0 


10000 


7500 


5200 


2.5 


1.5 


15 


75 


0.16 


1.0 


Q0842+3431 


2.100 


14.0 


5.00 


16000 


10000 


8000 


3.0 


1.5 


15 


78 


0.04 


1.0 


Q1235+1453 


2.835 


5.00 


10.0 


13000 


8000 


7000 


2.5 


1.5 


15 


76 


0.12 


1.0 


Q1239+0995 


1.970 


3.00 


5.00 


16000 


8000 


8000 


3.0 


1.0 


15 


82 


0.10 


1.0 


Q0019+0107 Ml 


2.075 


1.00 


15.0 


16000 


9000 


8000 


2.5 


1.5 


15 


90 


0.11 


1.0 


QOO 19+0 107 M2 


2.075 


0.00 


15.0^ 




9000 


8000 


2.5 


1.5 


15 


90 


0.11 


1.0 


Q004 1-4023 Ml 


2.450 


5.00 


10.0 


10000 


5000 


3000 


2.0 


1.5 


15 


90 


0.18 


1.0 


Q004 1-4023 M2 


2.450 


15.0 


5.0 


6000 


3000 


3000 


2.0 


0.5 


15 





0.07 


1.0 



" kp„, represents the total equatorial (i = 90°) optical depth integrated over frequencies since t,o,^ = 0. 



5. Discussion 

We showed that a simple two-component equatorial -i-polar wind 
model is able to reproduce a variety of BAL profiles, ranging 
from detached absorption troughs to P Cygni-type profiles. The 
solutions of the fits are not unique and several models with dif- 
ferent geometries and/or physical properties can equally repro- 
duce the observed spectra. In accordance with previous studies 
(e.g. Hamann et al. 119931 ). this demonstrates that a unique phys- 
ical characterization of the outflow cannot be derived from line 
profile fitting. 

While detailed information on the geometry of the outflows 
cannot be derived, we nevertheless reached some interesting 
conclusions. First, in some objects, it is necessary to include 
both the equatorial and the polar absorption regions. This is in- 
deed the case for objects like Q1413-H117, Q1333-H2840, but 
also Q1235+1453 or Q0842 -1-3431, where the polar component 
allows reproduction of the shallow absorption trough observed at 
higher velocities and where the lower velocity equatorial compo- 
nent is needed to reabsorb the emission from the polar wind. In 
other objects (the prototype in our sample being Q0019-I-0107), 
the polar component is not necessary and the profiles can be fit- 
ted with only a rapidly rotating equatorial wind seen nearly edge- 
on. In this particular case, there is a strong degeneracy between 
the model parameters, as also suggested by similar profiles hav- 
ing been computed in the framework of other wind models (e.g. 
Proga i2003l Proga & Kallman [2004] l. 

Interestingly, in many of the fitted spectra, the viewing an- 
gle to the wind axis is found to be high, suggesting that BAL 
quasars are essentially observed when the optically thick equa- 
torial wind blocks the direct view to the continuum source. This 
agrees with the high inclination generally inferred to account 
for the spectropolarimetric properties of BAL QSOs (Schmidt & 
Hines [19991 Ogle et al. 1999). In this case, the polar wind could 
play the role of the extended scattering region at the origin of 
the polarization. We also showed that, when the profile is consti- 
tuted of a single deep absorption trough and a quasi-symmetric 
emission peak (e.g. Q0041-4023), the line can be produced in 
a two-component wind seen at low inclination (see the bottom 
right panel of Fig. |5]). This result underlines that some BAL 
troughs can be observed as the result of a polar outflow in a 
two-component wind with a large covering factor As mentioned 
in the introduction, a growing number of radio observations 
of BAL QSOs provide evidence for such two-component wind 
models, while recent hydrodynamical simulations show that a 
stable two-component wind can result from the flow around a 



black hole and its accretion disk (e.g. Proga l2007b . This model 
appears challenging for the unification by orientation scheme of 
BAL and non-BAL quasars (e.g. Turnshek ll984al Hamann et 
al. 1 19938 because of there is an absorption component at every 
viewing angle (cf. panel E of Fig. |2|. However, that the polar 
component can be omitted in some cases (e.g. Q0019+0107) 
supports the existence of a class of objects in which only a frac- 
tion of the continuum source is covered by the BAL material. 
All in all, the observations suggest that BAL outflows can have 
a wide range of covering factors. 

In our study we decided to use a simple radially expanding 
wind model with an equatorial and a polar component. This type 
of wind can account for most characteristics of the resonance 
line profiles observed in the spectra of BAL QSOs. However, 
we were not able to fit the C iv line shape of the P Cygni-type 
quasar prototype PHL5200. This indicates that the model used 
does not include all the ingredients needed. Indeed, the model 
cannot reproduce the very sharp transition observed between the 
absorption and the emission components in the Civ profile of 
PHL5200 (Turnshek et al. il988i ). Such a sharp transition at zero 
velocity could in turn be produced in a wind launched from the 
disk itself. In that type of model, which exhibits large-scale prop- 
erties similar to those of the model considered in the present 
study, the wind is launched from the disk and then radially ac- 
celerated by the radiation pressure (Murray et al. II995I Proga 
& Kallman 120041 ). When observed nearly edge-on, it can pro- 
duce strong absorption at the center of the line as observed in 
PHL5200. Proga ( 2003 ) shows the capabilities of these winds to 
produce a sharp transition between the absorption and emission 
in the case of cataclysmic variable stars. Implementing that type 
of wind is beyond the scope of this paper and has been left for 
future work. 

Furthermore, in a majority of objects we found it necessary 
to use high values of the ratio of the rotational speed to the po- 
lar terminal speed of the wind Vrot/Vmax to adequately reproduce 
the asymmetry of the emission line profiles. Such a high value 
of the rotation accounts for the redshift and the faintness of the 
emission peak. These characteristics result from the continuum 
emitted on the red side of the line profile possibly being ab- 
sorbed by the optically thick material located between the source 
of continuum and the distant observer when rotation is present 
(see Fig. [3]). The rotational velocity at the base of the wind can 
reach several thousand km s"', which is consistent with the ro- 
tational velocities inferred for gas orbiting in the vicinity of a 
supermassive black hole (e.g. Murray et al. II995I Young et al. 
,20071) . Interestingly, we found that the rotational velocity must 
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Table 3. Blueshift of the C iv line relative to the systemic redshift 
Zsys for the three quasars from our sample with secured [O iii] line 
measurement. 



BAL QSO Name Ref. z„,, Avbitieshift 

Q0019+0107 nil Dietrich et al. 2009 5500 km s"' 

Q0226-1024 2.268 Mcintosh et al. 1999 9200 km s"' 

Q1413+1143 2.553 Hutsemekers et al. 2009. 2500 km s'' 



Shields et al. l2009l l. Ultimately, it might be necessary to consider 
these effects for full self-consistent modeling of quasar outflows. 
In the particular case of Q0019+0107, two narrow Civ absorp- 
tion lines are observed at -3000 and +3000 km s ' in the wind 
rest frame (Fig. 4), i.e. at -8000 and -2000 km s"' in the sys- 
temic quasar+host rest frame. These velocities suggest that they 
originate in a large-scale outflow in the host galaxy rather than 
in the BAL wind. Interestingly enough, the velocity difference 
between these narrow absorption lines is close to the velocity 
separation between the Lya and NV resonance doublets (~ 5900 
km s"'), suggesting a possible line locking effect (Korista et al. 
IT9931 Arav & Begelman H994| l. 

6. Conclusions 

In this study, we used a combination of a Monte Carlo radia- 
tive transfer code and a simple two-component polar -i-equatorial 
wind model in which photons are emitted from a central spher- 
ically symmetric source and resonantly scattered in the wind to 
reproduce typical Civ resonance line profiles selected from a 
homogeneous sample of BAL QSO spectra. 

Although the lack of uniqueness of the line profile fitting 
does not allow us to strongly constrain the geometry of the wind, 
we can summarize our main findings as follows 

1 . The diversity of BAL profiles produced by the adopted po- 
lar+equatorial model ranges from the typical P Cygni-type 
profiles to the detached absorption ones, reproducing those 
observed in a homogeneous sample of BAL QSOs. 

2. While in some cases the line profiles can be reproduced by 
a single equatorial wind, we find it necessary to use a two- 
component polar -Hequatorial wind in a majority of objects. 

3. The viewing angle to the wind is generally large (disk seen 
near edge-on); however in some cases, the line profiles can 
also be reproduced when assuming a pole-on view, in accor- 
dance with the results of recent radio surveys of BAL QSOs. 
In this context, it would be interesting to obtain good quality 
spectra of bona-fide polar BAL quasars and try to fit their 
line profiles by assuming the pole-on view. 

4. The equatorial wind is rotating, and the rotational velocity at 
the base of the wind can reach a significant fraction of the 
polar terminal speed. 

A possible way to break the degeneracy between the various 
parameter combinations of the two-component model that can 
reproduce the observed BAL profiles is to use gravitational mi- 
crolensing. Indeed, a microlens moving across the quasar inner 
regions can differentially magnify the different line-forming re- 
gions, inducing line profile variations from which the geometry 
of the outflow can in principle be retrieved (e.g. Hutsemekers 
1993] Hutse mekers et aL [T994l Lewis & Belle 11^981 Chelouche 
2003] l2005l l. Our code MCRT has been explicitly built to inte- 
grate these microlensing effects. The effect of microlensing on 
BAL profiles, their use for deriving the physical properties of 
the outflow, and application to a known lensed system will be 
presented in a second paper. 
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Fig. 4. Illustration of the best fit for eight out of the ten Civ line profile selected from the BAL quasar sample of Korista et al. 
( I1993I I. The observed spectrum is represented by a full line and the two-component model profile by a dashed line. Each spectrum 
is normalized to the continuum level and represented as a function of the velocity. 
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Fig. 5. Illustration of the non-uniqueness of the fit for two particular objects of our sample. The upper panels (labelled Ml) of this 
figure show that the C iv line profile of these two objects can be fitted in the framework of the two-component wind model. In the 
lower panels (labeled M2), we show that the lack of typical two-component wind signatures in the line profile allows alternative fits 
with different parameters (see text). As in Fig.|4]the observed spectrum is represented by a full line and the two-component model 
profile by a dashed. Each spectrum is normalized to the continuum level and represented on a velocity scale. 



